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Abstract. Local field effects on the rate of spontaneous emission and Lamb shift 
in a dense gas of atoms are discussed taking into account correlations of atomic 
center-of-mass coordinates. For this the exact retarded propagator in the medium is 
calculated in independent scattering approximation and employing a virtual-cavity 
model. The resulting changes of the atomic polarizability lead to modifications of 
the medium response which can be of the same order of magnitude but of opposite 
sign than those due to local field corrections of the dielectric function derived by 
Morice, Castin, and Dalibard [Phys.Rev.A 51, 3896 (1995)]. [] 



1 Introduction 

The experimental progress in cooling and trapping of atoms and the ob- 
servation of Bose-Einstein condensation in atomic vapors has lead to a 
growing theoretical interest in the interaction of light with dense atomic gases 
HOI'il'll' 111 t ne present paper it is analyzed how the spatial distribution of 
nearest neighbors in a dense gas, characterized by two-particle correlations, 
affects the interaction of an excited atom with the electromagnetic field vac- 
uum. In particular modifications of the rate of spontaneous emission and the 
Lamb shift are calculated and the resulting modifications of the dielectric 
function discussed. 

While the interaction of light with a dilute gas is well described in terms 
of macroscopic quantities with the well-known Maxwell-Bloch equations, this 
is no longer true when the gas becomes dense. Here two new types of ef- 
fects arise: When the density of atoms becomes large enough such that the 
resonant absorption length is less than the characteristic medium dimen- 
sion, re-absorption and multiple scattering of spontaneous photons need to 
be taken into account. Secondly with increasing density dipole-dipole inter- 
actions between nearest neighbors become important. Here the macroscopic 
picture of a homogeneous polarization breaks down and it is necessary to 
introduce local field corrections. 

The most famous local-field correction leads to the Lorentz-Lorenz (LL) 
relation between atomic polarizability a(tu) and dielectric function e(uj) 
The LL-correction removes the unphysical contact interaction of two atoms at 

1 This paper is dedicated to the memory of Dan Walls 
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the same position which arises in a continuum picture of a homogeneous po- 
larization and is independent of any specifics of the atoms. Recently Maurice, 
Castin, and Dalibard || have derived a generalization of the LL relation that 
takes into account center-of-mass correlations for finite distances at which 
point specific properties of the atomic gas enter. They showed in particular 
that the tendency of bosonic atoms to bunch, when the critical tempera- 
ture of condensation is approached, leads to a measurable change of the 
complex refractive index. Even more pronounced effects such as a dramatic 
line-narrowing were recently predicted by Ruostekoski and Javanainen for a 
Fermi gas in the case ideal case 1 10 or in the presence of a BEC transition 



0- 

On the other hand, it is known since the early work of Purcell |12j , that the 
microscopic environment of an excited atom can also change its interaction 
with the field vacuum. Thus local field effects should lead to a modification 
of the atomic polarizability itself, in particular the spontaneous emission 
rate and the Lamb shift. Different macroscopic models for Lorentz-Lorenz- 
type corrections of the spontaneous emission rate of an atom embedded in 
a lossless dielectric have been developed jl3 14 l^ , |il||l7| and experimentally 
tested @@. In the presence of losses, as is the case in dense gases of the 
same kind of atoms, macroscopic [§0|,^l| and microscopic approaches [||| 
have indicated however that local-field effects due to nearest neighbors may 
be equally important as the LL correction of the contact interaction. In the 
present paper I analyze these effects using a Greens function approach. 

According to Fermi's golden rule the rate of spontaneous emission _T and 
the Lamb shift A are given by the (regularized) exact retarded propagator 
D of the electric field at the position t*o of the probe atom p2 23 24[. 



r = d ■ Re 



A = d ■ Im 



D(r ,r ;wo) 
D(r ,r ; w ) 



■d, 
■ d. 



(1) 
(2) 



d is the dipole vector and u>q the true transition frequency of the atom. From 
the known propagator in free space and after regularisation one finds for an 
isolated atom 



r = r n = 



A = 0. 



(3) 



The exact retarded propagator in a medium can formally be obtained from a 
scattering series. This series is here calculated neglecting multiple scattering 
of photons by the same atoms and taking into account only two-particle 
correlations of center-of-mass coordinates. It is shown that the alterations 
of the atomic polarizability a due to local-field corrections of spontaneous 
emission and Lamb shift lead to modifications of e which are of the same 
order of magnitude as those found by Maurice, Castin and Dalibard || and 
Ruostekoski and Javanainen pO| , pT| . 
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2 Scattering series for the retarded propagator 

A dense medium affects the interaction of an excited probe atom with the 
surrounding electromagnetic vacuum by multiple scattering of virtual pho- 
tons emitted and re-absorbed by that atom. These scattering processes can 
formally be described by the exact retarded Greens-function (GF) 

D(ri > r 2 ;r)=«(r)<0|[i;(ri,t 1 ),B(r a) i 2 )]|0), (4) 

where r = ti — t% and E is the operator of the electric field interacting with 
all atoms. The free-space or vacuum GF in the frequency domain is given by 
25U26 



x o x 



1 



pi(fc+iO)x r 

G^(x,u) = -k 2 P (ikx) 1 + Q (ikx) +-S(x)l, (5) 

47ra L x J 3 

where D° = i?iG /eo. Here k = lo/c, x = \x\ and 

P(z) = l-- + \, Q{z) = -l + --\. (6) 

z z z z z* 

G (x, oj) diverges as x — > which is related to the large-q behavior in recip- 
rocal space. This will lead to corresponding divergences in the exact propa- 
gator G(x, oj) which can however be removed by introducing a regularisation 
G(q,uj) — > G(<7, uj) f(A, q) with a wave-number cut-off A. For the purpose 
of the present paper I will assume that the large-q behavior of the Greens 
function is properly regularized and ignore all contributions containing the 
regularisation parameter A. The subject of regularisation will be discussed in 
more detail at a different place. 

The net effect of all possible multiple scattering events can be described 
by a scattering series for the exact propagator. We here assume not too large 
densities, such that dependent scattering can be neglected. I.e. there can be 
as many scattering events as possible but never twice from the same atom. 
In this so-called independent scattering approximation (ISA), the scattering 
series can be expressed in the form 

G(0,u/) = G°(0,cj) - ^G°(x 0i ,u/) • Bi(w) • G ^,^) + 

+ G°(x 0l ,uj)-!L l (uj)-G (x l] ,uj)-a :j (uj)-G (x j() ,uj) + ---, (7) 

where a.j(u>) is the polarizability tensor of the jth atom of the host material, 
the summation is over all atomic positions, and cc,y = — Ta. It should be 
noted that the polarizability of the excited probe atom does not enter the 
scattering series. However, the probe atom can affect the spatial distribution 
of the surrounding scatterers and it is necessary to keep track of its presence. 

We now assume a homogeneous medium of density g with randomly ori- 
ented two-level atoms such that a,i(u>) — a(uj) 1. In this case we may replace 
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the sums over atomic positions by integrals. For this we introduce normalized 
joint probabilities ^(^Oi r i), P3( r o, r i, r 2) etc. to find one atom at position 
ri and the probe atom at t*o; to find two atoms at positions n and r 2 and 
the probe atom at ro etc, with pi(r) = 1. This leads to 

G(0) = G°(0) -got J d 3 r, p 2 (r (h r i )G (x <H ) ■ G°(x i0 ) + 

+g 2 a 2 f fd 3 ri d 3 rj p 3 (r Q , r t , r 3 ) G a (x Qt ) ■ G°(x tJ ) ■ G°(x j0 ) + ■ ■ ■ (8) 



where we have suppressed the frequency argument for notational simplicity. 

For a dilute gas the positions of the atoms can be treated as independent 
and one can factorize all particle correlations, which amounts to p m = 1. 
The scattering series can then easily be solved. The poles of G(q,uj) = 
Jd 3 x e lq x G(x,oj) determine the dielectric function for which one finds the 
well-known dilute-medium result 

e(k) = f| = 1 + ga(w). (w - kc) (9) 

Properly regularizing G(q, uj) and transforming the result back to coordinate 
space eventually yields the spontaneous emission rate and the Lamb shift 
relative to the vacuum 

r = r Re [Vi] , A = -j- Im . (10) 

Here e is the dielectric function of the gas at the true transition frequency. 
For an atom embedded in a dielectric host with real dielectric function the 
result for r is identical to that obtained by Nienhuis and Alkemande based 
on a density-of-states argument 

3 Local field effects and center-of-mass correlations 

In order to obtain a non-perturbative result for the retarded propagator in a 
dense gas, one has to sum all contributions of the scattering series (||) with- 
out factorizing the center-of-mass correlations. Apart from some special cases 



such as hard-sphere scatterers 27 it is not possible to bring the scattering 
series in an exact closed form and approximations are needed. The first ap- 
proximation I use here is to take into account only two-particle correlations 
by using a Kirkwood-type factorization p8| of higher-order contributions 

P3(l, 2, 3)= pa(l, 2)p 2 (2,3)p 2 (l, 3), etc. (11) 

Furthermore it is assumed that only correlations between successive scatterers 
matter, which is however correct up to second order in the density. With this 
the scattering series (Js[) can be represented in the form 

G(0) = G°(0) - gaH a (x 0l ) ■ G a (x lQ ) + 

+g 2 a 2 U°(x Ql ) ■ U°(x tJ ) ■ W°(x j0 ) - (12) 
~g 3 a 3 H°(x Qi ) ■ H°(x tJ ) ■ W°{x jk ) ■ H°(x kQ ) + ■■■, 
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where the spatial integration has been suppressed and 

H°(r 1) r 2 ,a;)=p 2 (r 1 ,r 2 )G (ri,r2,w) (13) 

is the retarded propagator modified by the two-particle correlation p 2 . Note 
that the first order term contains only a single function H . 

Since two-particle correlations between atoms of the host material can be 
different from correlations between the (excited) probe atom and a host atom, 
we will distinguish these two in the following. This also includes the case of 
an atomic impurity in an environment of a different species. The scattering 
series can then be written in the form 

G(0, w) = G°(0, w) - ga /d 3 nH°(r , n, w) • G°(n, r„, w) + 

/d 3 n d 3 r 2 H°(r , n, u) ■ T< 2 ) (n, r 2 , w) • H»(r 2 ,r , W ), (14) 



where T^ 2 ) (n, r 2 , w) is the part of the scattering matrix that contains at 
least two scattering processes in the host material. is the free propagator 
modified by the correlation of the excited probe atom with an atom of the 
background gas. Although in ISA the T-matrix does not contain scattering 
events from the probe atom, it would in general still depend on its presence 
through the position correlations. With the earlier assumption that only cor- 
relations between successive scatterers matter, this dependence is lost. I.e. 
we treat the scattering of photons in the gas as if the place of the probe atom 
would be filled with the host material, which is equivalent to the virtual cav- 
ity model of Knoester and Mukamel |l4j]. From the discussion of impurities 
in cubic dielectric host materials by deVries and Lagendijk jlTj it is however 
expected that this approximation does not affect the results in leading order 
of the density. Finally we will restrict ourselves to the case of a homogeneous 
and isotropic gas, such that the two-particle correlation depends only on the 
distance between the atoms p 2 {r\,r 2 ) = p 2 {x\ 2 ) where x\ 2 = |ri — r 2 |. In 
this case the scattering matrix obeys a simple Dyson equation in reciprocal 
space 

T(q,w) = Qa(u)U g {q,u) Q a{u) - ga(u)H.° g (q,Lu) ■ T(q,u) (15) 

H" is the free propagator modified by the two-particle correlations of the 
host material. 

It is convenient at this point to introduce the irreducible correlation h 2 
according to p 2 = 1 + h 2 ■ One then has 

- - i k 2 it 9 + m*») i + g fe M > q -f . as) 



with 



/f (q,u) = -k 2 f^dxxe^ hti(x)\j (qx)P(ikx) + j -l^-Q(ikx)} , (17) 
J I qx J 

p OO 

f%(q,uj) = k 2 dxxe ikx h 2 l {x)j 2 (qx)Q(ikx). (18) 
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j n (z) are spherical Bessel functions and P and Q have been defined in eq. (g) . 
In the first term of eq.(|l6|) we have made use of p% (x) 5(x) — which cor- 
responds to the LL correction of the contact interaction. The solution of eq. 
( p"5| ) is now easy to obtain. The poles go of T which determine the dielectric 
function follow from the equation 



% 



QOt{tjj) 



/f(go,w)(«o 



iO = 0. 



(19) 



Since in lowest order of ga one has q 2 , 
ff(k,u>) which yields 



e(k) = 



4 
k 2 



1 



ga{u>) 



1 — ga(uj)/3 + ga(<jj)fi{k,u>) 



k 2 + iO, we may replace f((qo, oj) by 



(20) 



This result is identical to that of Maurice, Castin and Dalibard Q in ISA. It is 
interesting to note that although the free GF H° contained also longitudinal 
components (proportional to Jf q°q/g 2 ), they exactly cancel in the expression 
for the dielectric function. 

To obtain the spontaneous emission rate and Lamb shift, we consider only 
the leading order corrections in the density where we can replace /f 2 (<3S w ) by 
/(* 2 (fc,w) and introduce a regularisation of the large-q behavior. This yields 
after some algebra for the orientation averaged retarded GF: 



G(x = 0,w) 



Gtt 



1 + ga 



7 



ft 



-Q or ff 

24 3 Jl 



ft 



e 2 



If. 



2/r/f 



(21) 



One recognizes that there is again no contribution from the longitudinal 
terms / 2 up to second order in g. Furthermore the two-particle correlation 
ff between ground-state atoms enters only in second order of the density, 
while there is a first-order contribution from ff. This is physically intuitive 
since correlations between ground-state atoms enter only after two scattering 
events, while correlations involving the probe atom are already important in 
first order. From the above results one finds in leading order of the density 



r = r a 



1 + \&ol - Q (a'ff - a" ft") 
o 



2 



ga 



II fel i / rel 

g{a / x +af x 



(22) 
(23) 



where we have introduced the real and imaginary parts of ff = ff' + iff . 
Eqs.(22) and ( p3| ) are the main result of the present paper. The first-order 
corrections to spontaneous emission rate and Lamb shift, which are inde- 
pendent on the /'s are the Lorentz-Lorenz local field corrections derived in 
[E2|. 
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In order to illustrate the implications of eqs.(g2|) and ( p3[ ) to the medium 
response, I will discuss in the following section the dielectric function of a 
dense gas of two-level atoms using some simple model functions for the center- 
of-mass correlations. 



4 Modifications of medium response 



We have shown that local field corrections change not only the relation be- 
tween atomic polarizability and dielectric function of the medium but also 
the atomic polarizability itself (with respect to the thin- medium case). To 
illustrate the net effect of both corrections to the dielectric function, I now 
consider a classical, homogeneous gas of radiatively broadened (cold) two- 
level atoms with randomly oriented dipole vectors. The dimensionless atomic 
polarizability of such atoms is isotropic and has in free space the strength 



a (d) = ak^ 



(in 



(24) 



where kg is the resonance wavenumber, S = (uj a b ~ ^)/lab is the normalized 
detuning from the true resonance and 7 a f, = i~b/2 is the free-space dipole 
decay rate. With this we find in lowest order of the dimensionless density 

q = g/kl 



r = r Q 



2 



7ir g — 6tt g /f 



(25) 



Substituting these expressions into the dielectric function, eq.(pfj|), one even- 
tually finds for the shift relative to the dilute-medium resonance and the 
effective linewidth up to first order in g: 



7cff = T 



{o 
2 

oft 2 



1 - 67rp/f" + 6irgfi" 



—2irg + 6irg ff' + 7tt g — 6ir g ff' 



(26) 
(27) 



The second term in the expression for the linewidth is due to local field 
corrections of the dielectric function and the third one due to the changed 
spontaneous emission rate. It is interesting to note that both contributions 
are of the same order of magnitude but differ in sign. Thus local field effects 
to the vacuum interaction may compensate the line-narrowing / broadening 
effects resulting from local field corrections of the dielectric function in lowest 
order of the density. In the expression for the line-shift, eq.(|27]), one recognizes 
the familiar Lorentz-Lorcnz shift —2np. The second term emerges again from 
local field corrections of the dielectric function and the two last terms are 
due to modifications of the Lamb shift. 

To illustrate the effect of the center-of-mass correlations let us consider 
a gas with repulsive interaction such that the two-particle correlation /12 is 
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close to —1 over an effective correlation distance z and then approaches zero. 
As simple model functions we use a Gaussian and a hyper-Gaussian 

h { 2 a) (x) = - exp [-x 2 /z 2 ] and (x) = - exp [-x 8 /z 8 ] (28) 

where the distance x, — x/X and the correlation length z — z/X are normalized 
to the resonance wavelength A. I have plotted in the following figures the 
real and imaginary parts of /i for both correlations. It is worth noting that 
for values of the correlation length larger than the resonance wavelength, 
/' — > 7/12, while /" becomes a linear function of z. 



0.1 




0.5 1 1.5 2 2.5 3 2 4 6 8 10 



z z 

Fig. 1. Real (/') and imaginary parts (/") of / as function of normalized correlation 
length z — z/\ for Gaussian and hyper-Gaussian correlation functions and 



If the correlation length z scales with the density according to z ~ pT 1 !^ 
the above behavior can be associated to the density dependence of the correc- 
tion terms. In the following plots the real and imaginary parts of z~ 3 fi <~ gf\ 
are shown as function of z ~ 3 ~ g. One recognizes that gf" scales for small g as 




g 2 / 3 j while gf[ is approximately linear in the density. Thus in the low-density 
limit linewidth changes will dominate line shifts. 
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5 summary 

In the present paper local field effects on spontaneous emission and Lamb 
shift in a dense atomic gas have been discussed taking into account two- 
particle center-of-mass correlations. It has been shown that the corresponding 
changes of the atomic polarizability can lead to modifications of the medium 
response, which are of the same order as those resulting from direct cor- 
rections of the dielectric function found by Maurice, Castin and Dalibard 
|1 and Ruostekoski and Javanainen [fiof . They are however of opposite sign 
and thus may compensate the leading order corrections in the density. The 
present approach is based on an independent scattering approximation and 
a virtual cavity assumption. It is thus only applicable for densities which are 
sufficiently smaller than the cubic wavenumber, p« 1. Also the existence of 
a distinguished probe atom has been assumed, which is valid only for classi- 
cal gases. Nevertheless the results indicate that local field corrections to the 
atomic polarizability may change or even reverse the predicted line-shifts and 
linewidth modifications found for Bose gases near the condensation temper- 
ature and for low-temperature Fermi gases. 
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